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Abstract 

We demonstrate with soluble models how to employ the effective Hamil- 
tonian approach of Lee and Suzuki to obtain all the exact eigenvalues of 
the full Hamiltonian. We propose a new iteration scheme to obtain the 
effective Hamiltonian and demonstrate its convergence properties. 



1 Introduction 



Much effort [1-5] has been made in the past few decades to calculate the shell- 
model effective interaction from a realistic nucleon-nucleon (NN) interaction, 
which reproduces properties of the deuteron and the NN scattering data. A 
completely satisfactory solution to this problem has not yet been found de- 
spite much success || 0. One of the major difficulties is that the Rayleigh- 
Schrodinger perturbation series for the full effective interaction in terms of the 
nuclear reaction matrix (G-matrix) (s) does not show convergence within the 
first few orders, and the degree of effort required to accurately calculate the 
higher-order terms grows rapidly In fact, Schucan and Weidenmiiller [fiof 
have shown that in cases when intruder states are present in the model space, 
the perturbation series must diverge. 

Another problem which remains largely unsolved is that shell-model results 
are found to be dependent on quantities which are, in principle, arbitrary. For 
example, the results are found to be dependent on the unperturbed single- 
particle Hamiltonian, Ho, and on cj, the starting energy of the G-matrix. The 
residual dependence of the results on these quantities is a measure of the lack 
of a fully converged effective Hamiltonian. Of course, one may select a self- 
consistent Hq to eliminate classes of higher-order diagrams and one may choose 
u> to give reasonable agreement between theoretical and experimental spectra. 
However, a fully converged effective Hamiltonian would have a strong measure 
of independence of both Hq and u>. 

Some ten years ago, Lee and Suzuki (ll[] proposed an approach to computing 
the effective interaction based on a similarity transformation. They suggested 
two iteration methods which are able to assure the independence of the results 
from the starting energy and, therefore, to bring a degree of convergence to 
the effective interaction. To employ their approach, however, it is necessary to 
evaluate the "Q-box" , an infinite series of irreducible, valence-linked diagrams 
[| U). Thus, the problem still remains to sum the series that defines the Q-box 
to all orders. 

The effective Hamiltonian is greatly simplified in a "no core" model space 
[ fl3| , thereby avoiding some of the convergence difficulties. Since there is no inert 
core, there are no hole lines and the Q-box reduces, at the two-particle effective- 
interaction level, to the bare G-matrix, which has only the ladder diagrams 
(to all orders), and can be calculated exactly |L4[ [Hi). The only corrections, 
in principle, are many-body forces, but first one must solve the problem for 
effective two-body forces accurately before meaningful statements can be made 
about these higher-body forces. In the no-core approach, all the particles in a 
nucleus are active, so for practical reasons one has to restrict oneself to light 
nuclei. 

Poppelier and Brussaard [jl6| have adopted a no-core model space, evaluated 
the G-matrix as the leading contribution to the Q-box, and have applied the 
Lee-Suzuki method to light nuclei. Their work provides some motivation for the 
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present endeavor, since they obtained encouraging results, even though they did 
not achieve ^-independence of the effective Hamiltonian. We are able to obtain 
w-independent effective Hamiltonians in our applications to soluble models. 

We will demonstrate two additional properties of the Lee-Suzuki approach: 
independence of the effective Hamiltonian from the choice of the single-particle 
Hamiltonian H Q and the ability to obtain all the exact eigenvalues of the full 
problem. These demonstrations will be given with soluble models. 

The organization of this paper is as follows: We first give a brief review of the 
Lee-Suzuki formalism and its application by Poppelier and Brussaard. We show 
the independence of the resulting effective Hamiltonian from Hq. We then sug- 
gest a different iteration equation for the effective interaction than the one used 
in Ref. |]l6| and show that when the truncated model space is one-dimensional, 
this iteration equation can be reduced to a simple, intuitively correct expression, 
stating that the total energy is equal to the unperturbed energy plus the inter- 
action energy. This leads to the demonstration that all the converged effective 
energies are exact eigenenergies of the system. We finally apply our procedure 
to some simple systems, investigate the convergence properties of the approach, 
and show that all the exact energies are obtained in the truncated model space. 



2 H e g: Method and Properties 

In this section, we give a brief review of the Lee-Suzuki formalism JO], [T^, [l6| 
and adapt it to our version of the nuclear no-core problem, point out the Ho 
independence and introduce an alternative iteration scheme. 
Let us consider the eigenvalue problem 

HV = EV, (I) 

where H is the full Hamiltonian and \& is the wave function of an A-particle 
system. The Hamiltonian H consists of the relative kinetic energy term T le i 
(i.e., the center of mass kinetic energy has been removed) |f7| and the potential 
term V, which we take to be a purely two-body operator: 

H = T Icl + V = H + (T rel + V - H Q ) = H + H u (2) 

where Ho is the arbitrary unperturbed Hamiltonian and Hi = H — Ho- 

We introduce the operator P which projects onto the model space and the op- 
erator Q which projects onto the excluded space, so that, 1 = P+Q, PQ=QP=0 
and [Ho,P] = [H o ,Q]=0. 

We then consider a similarity transformation: 

H = e~ s He s 

= PHP + PHQ + QHP + QHQ, (3) 
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where S has the property that S = QSP (hence, PSP = QSQ = PSQ = 0, 
S 2 = S 3 = ■ ■ ■ = 0, e s = 1 + S, e- s = 1 - S). For PHP to be a P-spacc 
effective interaction P c ff , we require 

QHP = Qe- s He s P = 0. (4) 

It then follows that if S exists, any P-space eigenvalue is also an eigenvalue of 
H and, therefore, of P. Equation (|j), which determines S, can be rewritten as 

QP/P + QHQS - SPHP - SPHiQS = 0. (5) 

After this equation is solved for S, the effective Hamiltonian is given by 

P off = PHP = PHP + PP/QS, (6) 

and the effective interaction by 

V cB = P c ff - PP P = PP/P + PP/QS. (7) 

It has been shown by Poppclier and Brussaard |l6| that Eq.(||), after a 
quantity ojS is added to both sides, can be rewritten as 

S = 7TWn QHlP ^rr^SPiHo-oj + V^P. (8) 

uj — QHQ uj — QHQ 

In order to solve the above equation iteratively for S, one introduces a quantity 
Z: 

Z = PH a P-uP + V cSl (9) 
and defines the generalized G-matrix: 

G{uo) = PHiP + PHiQ L—Q^P. (10) 

u> — QHQ 

This generalized G-matrix (referred to as the "Q-box" in Ref. |LlJ) is an A- 
body operator since P and Q project onto ^4-particle states. In applications to 
nuclei, one expects the Brueckner-G-matrix B to be a leading approximation 
to Eq.(10), and this is the approximation invoked in Ref. In our model, we 
retain the complete generalized G-matrix throughout. 
Next one rewrites Eq.(|J) as 

(1 + PP/Q L^ S )Z = PH P-uP + G(co). (11) 

Poppelier and Brussaard suggest an iteration equation for the solution of Z 
{Zt = 0,n > 2): 

z n = -j — ^ pr~^ ~ n i^T? ? P{H -uj+Go)P 

(12) 
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where Go=G{u>) and Gk is related to the /c-th derivative of G(u>): 



Gk{Ld) -H^U^- (13) 

If the iteration procedure converges, the effective interaction can be obtained 
from Eq.(g). 

One can now see that H e g = Z+ujP is independent of Hq since G(u) depends 
only linearly on Hq, leading to the cancellation of the Hq term in Z n , Eq.(|l2|). 
This manifest independence of Hq is one of the more appealing features of this 
method. 

In any practical application of the iteration procedure (|l2|), we can only 
include a finite number of derivatives Gk- Assuming Gk — for k > N, we can 
rewrite (12) for n > N as 

z n = "j ^ pri 7 7=^7 ~ 7^-^ ~ ~ P(i7 -^+G )P. 

1 — Cxi — Cx2^n-l — Cx3Z„_2Z„_i — • • • — Cxjv^ri-A r +l^n-AT+2 1 ' ' <&t»-l 

(14) 

It is easy to see that if the iteration is to converge, the converged value Z has 
to satisfy the following equation which can be obtained from Eq.(|lJ) by taking 
the limit n —> oo and identifying Z n=OQ = Z: 

z = x- Gi -g 2 z-gIz^---g n z^ p ^ - u + G ^ p - ™ 

We, therefore, obtain an alternative iteration equation: 
Zn = ——-— pr- —^[Hq-u + G^P. (16) 

-L — Cxi — tx2^ri-l ^ Cx3^ n _ 1 — • ■ ■ — LriV/5 n _ 1 

We should point out that the above iteration equation can also be applied to 
a case for which the k-th derivative of the G-matrix does not vanish (but is 
sufficiently small so the iteration converges) for any large k. 



3 Test Case: One-Dimensional Model Space 

We now show that for the case in which the model space consists of only one 
unperturbed basis state (so all the matrices in the model space become num- 
bers), Eq.(15) can be reduced to a particularly simple form which is intuitively 



appealing. In fact, we have for N — > oo, that 

Go + G X Z + G 2 Z 2 + ■■■ + G N Z N G{uj + Z). (17) 



So Eq.(|15[) becomes (Eq is the eigenenergy of the unperturbed Hamiltonian 
PH Q P) 

z = i + tGH-U + s)]/* [E °~" + G(w)] (18) 
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which can be easily simplified to yield 



uj + Z = E + G(u + Z). (19) 

This is an equation which a converged value of Z must satisfy for any starting 
energy ui. It is evident from Eq.([l9]) that the effective interaction, given by 
T4ff = u> + Z — Eo, is independent of ui. In other words, Eq. (|Ts|) shows that for 
any chosen u>, Z adjusts so that u> + Z is independent of that u. Remembering 
that ileff — o-> + Z, we note that Eq.(|l9|) is simply stating that the total energy 
E=H e g is the unperturbed energy Eq plus the interaction energy V c s, or, 

E = E Q + G(E), (20) 

after we identify V c g — G(E). Note that E again is independent of H (or E ) 
since G(E) has a term, -E , cancelling the explicit E on the right-hand side of 
Eq.©. 

Let us consider a simple example for which G{uj) has a quadratic dependence 
on lu: 

G(u) = (g - E ) + Bl u + -g 2 uj 2 , (21) 

where go, gi, and g 2 are constants independent of lo. We have for this choice of 
G(u) that 

G = G(u>), Gi=g 1 +g 2 u>, G 2 = 52/2, 



and Gfc = for k > 2. Eq.(15) for Z then becomes (note cancellation of Eo) 

Z= l-Gi-G 2 Z [Eo ~ U> + G( -° j)] - (22) 
The solutions (denoted by Z + and ZJ) to this equation are 
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Hence, the resulting energies, given by E = H e s = H n + V e g = Z± + ui, are 
clearly independent of the starting energy ui and of Ho- It is easy to verify that 
these energies E also satisfy Eq. (p]|) . 

We claim, as will be shown explicitly in the next section for a couple of 
simple cases, that with the standard definition Eq. (|Io|) for the generalized G- 
matrix, Eq. (|2(]) for the energy E in the model space is equivalent to the secular 
equation for the eigenenergies of the system in the entire space: 

det(H -EI)= 0, (24) 

where H is the full Hamiltonian in the matrix form and / is the unit matrix. This 
will then be sufficient to conclude that we may obtain all the exact eigenvalues 
by these procedures. 
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4 Application to Exactly Soluble Models 



In this section, we apply Eq. ( |20| ) , an outcome of the Lee-Suzuki iteration Eq. ( |l2| ) 
in the case of one-dimensional model space, to two exactly soluble models. We 
first consider a simple system with only a pair of unperturbed basis states, |1) 
and 1 2). Our model space is the space spanned by the state vector |1) and the 
excluded space is spanned by |2): P = |1)(1|, Q = |2)(2|. The full Hamiltonian 
H is then a 2 x 2 matrix and the unperturbed Hamiltonian Hq is any 2x2 
diagonal matrix: 



H 



h u 
h-21 



h\2 
h 2 2 



H 



a 
b 



The corresponding generalized G-matrix, defined by Eq.(|lO|), is then 

1 



G(u>) = (fin - a) + hi2 

So the equation for the effective energy, Eq.(| 

hi2h 2 \ 



, -h 2 i- 

becomes 



E = h u 



E-h 



•12 



or det(# - EI) = 0, 



(25) 



(26) 



(27) 



where / is the 2x2 unit matrix. 

Similarly for a system with three unperturbed basis states |1), |2), |3). We 
have P = |1)(1|, Q = |2)(2| + |3><3[ . The full Hamiltonian and the unperturbed 
Hamiltonian are written as 



(28) 





hu 


hi2 


hi3 




a 





" 


H = 




h,22 


fl23 


, Ho — 





b 







_ h 3 i 


h 3 2 


_ 










c 



The generalized G-matrix is 
G{uj) = (h n - 



*,J=2,3 



QHQ 



hn 



(29) 



with 



QHQ 



10 - h 2 2 

-h 32 



-h 2 3 
u - h 33 

1 



uj - h 33 
h 32 



h 2 3 
uj - h 2 2 



(30) 



or 



G{uj) = (hu-a) + 



(uj - h 2 2)(u - h 33 ) - h 23 h 32 
hi 2 (uj - h 33 )h 2 i + hi 2 h 23 h 3 i + ft.13/132/121 + hi 3 (u> - h 2 2)h 3 i 



(uj - h 2 2)(uJ - h 33 ) - h 23 h 32 



(31) 
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So Eq.(]2C|) again becomes 

det(H-EI) = 0, 



(32) 



and independence of Hq and u) is again guaranteed. 

Thus, we have seen by explicit construction, for a one-dimensional model 
space, that solving for the exact effective Hamiltonian gives rise to the entire 
spectrum by finding all solutions to Eq. (^0|) . Although we have not proven this 
property for the general finite dimensional model space problem, we nevertheless 
believe it to be true. 



5 Convergence Properties 



In order to examine the different convergence properties of the iteration equa- 
tions ( |l2| ) and (|l^), we arbitrarily make the following choices for H and Hq: 



H 



1 -0.5 -0.25 
-0.5 3 -V2 
-0.25 — v/2 4 



10 
2 
3 



(33) 



The generalized G- matrix, given by Eq.([3l|), becomes 

G{uj) = + — — - 



(34) 



with a = (9 + 4\/2)/48 and = (3 - 2a/2)/24. It has two simple poles at u=2 
and ui=5. Note that the poles of G(u>) are controlled by H and are independent 
of Hq. The exact eigenvalues E for this problem, which satisfy Eq.(|20|), or 



E=1 + G(E), 



(35) 



can be solved diagrammatically, as shown in Figure 1. In this figure, we plot 
the function w and the function 1 + G(uj). The solutions E are given by the 
intersections of the straight line u> and the curve segments for 1 + G{ui). 

In the same figure, we also show the converged value for the effective Hamil- 
tonian H c ff = Eq + V e {£ using the iteration equation Eq.(12) for a wide range 
of starting energies u> (dotted line). It is seen from the figure that all three 
solutions to Eq. (po|) can be achieved through the iteration procedure ( |l2| ) and 
that the solution which is reached depends on the starting energy that is chosen. 
When the starting energy u> is smaller than 2 (the position of the first pole) , the 
converged energy is identical to the lowest eigenenergy. When lu is in between 
the two poles, 2 < u> < 5, the first excited state energy is obtained and when u) 
is beyond the second pole, ui > 5, the second excited state energy is obtained. 
Note that the exact eigenvalues do not coincide with the poles of the generalized 
G-matrix. 



8 



The experience with this simple 3x3 example suggests that, in more realistic 
applications, one must sweep uo through a wide range of values to ensure that 
one obtains the ground-state eigenvalue. 

What is not evident from the figure is how many iterations are needed to 
achieve such a convergence. We show that in Table I. Our criterion for con- 
vergence is that the difference in the effective interaction between the final two 
iterations is smaller than 10 -7 , i.e., 

\V cS {N)-V cS {N ~\)\ < 1(T 7 , (36) 

where N is the number of iterations listed in Table I. ft should be pointed out 
that when Eq.(^|) is used, the number of derivatives (Gk) to be summed over 
in the denominator grows with the iteration number n ( i.e. /c max — fi — f ) but 
when the iteration Eq . (fl6|) is used, the number of derivatives is fixed. For the 
current application, we include Gk for k—l to fc max with fc max = 200. 

One sees from Table I that for a very wide range of the starting energy, 
it takes no more than 15 iterations to reach the converged value if Eq.(|16[) 



is used. When Eq.(12) is used, the number of iterations needed to achieve 
convergence varies with the starting energy and is generally larger than the 
number of iterations for Eq.([l6|). However, as also indicated in the table, when 



Eq.(16) is used, it is quite often the case that there is no convergence for the 
starting energy beyond the first pole u> = 2. 

It is also obvious from Table I that with the iteration Eq. (|l2|) , the number 
of iterations (hence, the number of derivatives of the Q-box) needed for the 
convergence specified by Eq.([56|) is small when the starting energy is close to 
the self-consistent value (i.e., the converged energy H e $). It grows rapidly as 
the starting energy is shifted away from the self-consistent value. 



6 Conclusions 

Through analysis and applications with soluble models, we have explicitly shown 
that the Lee-Suzuki approach to the effective interaction leads to a final result 
for the energy which is independent of the starting energy u> and independent 
of the choice of the single-particle Hamiltonian H . In addition, with soluble 
models, we show the method has the property that it can converge to any 
eigensolution of the full Hamiltonian. We have demonstrated this for 2x2 and 
3x3 matrix models; however, we believe it to be true for the general case of 
an N x N matrix. For a 3 x 3 matrix example, we explicitly indicate when the 
numerical result for a one-dimensional model space converges to the ground- 
state energy and when it converges to each of the excited-state energies. The 
state to which the method converges depends upon the value of the starting 
energy lu relative to the poles in G. We have also developed a modified iteration 
scheme which converges more quickly than the Lee-Suzuki approach but only 
in the case of the ground-state energy. These approaches hold great promise for 
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constructing accurate and meaningful effective-interaction matrix elements for 
shell-model calculations. Large scale calculations with the Brueckner G-matrix 
computed from different realistic free N-N potentials are in progress. 
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Table I. The convergence properties of the iteration procedures (|T^) and 
16| ) in a simple case for which the G-matrix is one-dimensional and of the form 
G(lo) given by Eq.(|34|). In the table, to is the starting energy which we vary. N 
is the smallest number of iterations needed to achieve the convergence defined 
by the inequality (|3^) . H e & = Z + to is the converged value for the effective 
Hamiltonian. For the iteration procedure ( |l6| ) we include the derivatives Gk 
to 200-th order (i.e., k = 1, 2, • • • , 200). The symbol "-" in the table means 
no convergence within 250 iterations. All three converged values for H e g are 
identical to the exact eigenenergies of the system. 



iV, 



Eg. (12) 



^Eq.(16) Hctt 



-20.0 
-10.0 

-5.0 
0.0 

0.75 
1.2 

2.25 
2.5 
3.0 
4.0 
5.1 
6.0 

10.0 

20.0 



177 
100 
59 
15 
4 
19 
5 
11 
18 
43 
8 
19 
51 
127 



12 
11 
8 
4 
9 
5 
54 



0.75338 
0.75337 
0.75337 
0.75337 
0.75337 
0.75337 
2.24479 
2.24479 
2.24479 
2.24479 
5.00183 
5.00183 
5.00183 
5.00183 



Figure Caption 

Figure 1. The converged effective Hamiltonian H c s = Z + oj, with G(lo) given 
by Eq.(f34|), as a function of the starting energy u> (dotted line) using Eq.(p"2|). 
The exact solutions for H c g are given by the intersections of the solid curves [lo 
and E Q + G(lj)], see Eq.©. 
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